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Abstract 

We study a question with connections to linear algebra, real algebraic 
geometry, combinatorics, and complex analysis. Let p(x, y) be a polyno- 
mial of degree d with N positive coefficients and no negative coefficients, 
such that p — 1 when x + y = 1. A sharp estimate d < 2N — 3 is known. 
In this paper we study the p for which equality holds. We prove some 
new results about the form of these "sharp" polynomials. Using these 
new results and using two independent computational methods we give a 
complete classification of these polynomials up to d — 17. The question 
is motivated by the problem of classification of CR maps between spheres 
in different dimensions. 

1 Introduction 

In this paper we answer by computational methods certain difficult questions 
about the set of bivariate polynomials with nonncgativc coefficients constant 
on a line. Following the notation of |DL09al IDLP07j . let H(2,d) denote the 
set of polynomials p(x, y) of degree d with nonnegative coefficients such that 
p(x, y) = I whenever x + y = 1. 

The condition that the coefficients are nonegative is motivated by a question 
in CR geometry, which we describe in § O Without this condition, the affine 
space I of all polynomials of degree d or less such that p(x,y) = 1 whenever 
x + y = 1 is easy to describe. It is in one to one correspondence with the 
vector space of polynomials of degree d — I or less. I.e. if q is of degree d — 1 
or less then we let p(x, y) — q(x, y){x + y — 1) + 1 el. When the coefficients 
of the polynomials are considered as variables, 1 is the solution set of a certain 
nonhomogeneous linear system. 

On the other hand, H(2, d) is a convex subset of I (with nonempty interior). 
It is the intersection of the positive cone, X, and the open set of polynomials 
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of degree exactly d. The questions we address are difficult because we need to 
consider the geometry of the boundary of %{2, d). When we try to answer these 
questions computationally, we note that the complexity grows very fast. We will 
describe two methods that allow one to effectively answer questions about those 
polynomials satisfying a certain extremal property. We prove certain theoretical 
statements of independent interest about these extremal polynomials, which 
allow us to reduce the complexity of computation. 

The condition that a polynomial is constant on a line leads to a system of 
linear equations, and hence to a problem in linear algebra. As we require the 
coefficients to be positive, we also get a natural formulation of our classifica- 
tion question below as a linear programming problem. Imposing an extremal 
condition we will describe makes it a mixed-integer programming problem. 

An important question from the point of view of CR geometry concerns the 
number of distinct monomials, N = N(j>), for p g H(2, d). N is the minimal em- 
bedding dimension of the associated CR map of spheres, see § [5] The following 
bound was proved in [DKR03 : 



are in fact polynomials in H(2,d) with d = 2N(fd) — 3, see |D'A 88 . Thus the 
inequality (fTJ) is sharp. 

The polynomials fd have many other interesting properties. In particular, 
they are group invariant. The CR maps that arise from fd are one of the only 
two possible classes of group invariant maps of balls. See § [§] for information 
and the references within. The polynomials fd are also related to Chebychev 
polynomials, arise in denesting radicals [OslOll, and have connections to number 
theory |D'A04| IDS05] . combinatorics |LWW04j and other fields. For example, 
fd(x, y) — x d + y d (mod d) if and only if d is an odd prime |D'A04j . 

If p €E 1-1(2, d) minimizes N(p) for a fixed degree d, then we call p a sharp 
polynomial. It can be proved that for a fixed degree only finitely many sharp 
polynomials exist. We wish to ask: Is fd the unique sharp polynomial in H(2, d) 
up to swapping of variables? If not, what are all the sharp polynomials for a 
given degree? 

The programs used for the computations were Mathematica 6 |Wol07j . and 
Genius 1.0.2 [Leb07aj . Some code was also written in native C for speed using 
the GMP [GMP07 library. We used two different approaches to the problem 
to get independent verification of the result and minimize the effects of possible 
bugs in the underlying code. The computer code used is available at the url: 
http : //www. jirka. org/LL08- archive . zip. 

Let us mention some relevant previous work. With the aid of a computer, 
Wono |Won93j classified all polynomials in H(2,d) with 5 terms or fewer, and 



d < 2N - 3. 
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Furthermore, for all odd d the functions 
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hence all sharp polynomials up to degree 7. The second author |Lic91| found 
new sharp polynomials in degree 11. D'Angelo and the first author |DL09aj 
have shown that there are infinitely many d for which there exist other sharp 
polynomials. 

The organization of the paper follows. In § [5] we state our main results 
including the list of all sharp polynomials in odd degree up to d = 17. In § [3] 
we prove new results about the form of sharp polynomials, which are useful in 
reducing the computation time. These results are of also independent interest. 
In §2] we describe treating coefficients of polynomials constant on x + y = 1 as a 
linear problem and prove several related results. In §[5] we describe the method 
of finding sharp polynomials by computing the nullspace of certain matrices. In 
§ [6] we describe the mixed- integer programming method we used to find sharp 
polynomials. In §[7] and §[8] we describe degrees in which uniqueness definitely 
fails by describing a construction of new sharp polynomials. This problem has a 
long and complicated history and motivation. In § [S] we discuss the motivation 
from CR geometry and complex analysis. 

We thank John D'Angelo for suggesting this research and for helpful dis- 
cussions along the way. We would also like to thank the referees for this and 
a previous draft for many useful comments and suggestions. The first author 
would like to acknowledge MSRI and AIM for holding workshops that focused 
on these and related questions. 

2 Main results 

In this section we state and discuss our main results. When p £ H(2,d) is such 
that the number of terms N(p) is minimal in H(2, d), we say that p is sharp. We 
say that uniqueness holds for degree d if there is exactly one sharp polynomial 
in H(2,d) up to interchanging the variables x and y. Otherwise, we say that 
uniqueness fails for d. Our main result is the following theorem. 

Theorem 2.1. Up to swapping of variables, Table\l\lists all sharp polynomials 
of odd degree for d < 17. In particular, uniqueness holds for d — 1, 3, 5, 9, 17. 

In the table, the first result listed is the group invariant one. To make the 
computations feasible, we prove new results about the form of sharp polynomials 
in § [3J In particular, in Lemma |3. II we prove that the degree d terms must be 
x d _|_ yd ^ an( j ^hese must be the only pure terms. We also prove that certain 
degree d— 1 monomials cannot arise, and that at least one degree d— 1 monomial 
is present. 

For degree d = 19 it is currently computationally infeasible to run the tests. 
Uniqueness does not hold in degree 19 by explicit construction (see § [7] and 
§ [5]) as 19 = 3 (mod 4) and 19 = 1 (mod 6). Hence, there are at least two 
inequivalent sharp polynomials in H(2,d) apart from the group invariant fd- 
When d = 21 it is unknown if uniqueness holds, although the method described 
in §[5] does not produce any new sharp polynomials for this degree. 
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To construct sharp polynomials of even degrees, we take two sharp polyno- 
mials of odd degree p and q, and write p(x,y) — x d + y d + po(x,y) utilizing 
Lemma |3~T1 We let 

f(x,y) :=x d + y d q(x,y)+p (x,y), (3) 

and we note that f(x, y) is sharp and deg / = degp + deg q. 

We have a computer assisted proof of the following theorem. Note that in the 
even degree case we can no longer make simplifying assumptions about terms 
of degree d or d — 1. Therefore, the computation is more expensive. The only 
assumption we can make is to use Lemma 13.61 to force exactly 1 pure term in 
each variable. 

Theorem 2.2. For 2 < d < 12, d even, the procedure in §3§ generates all sharp 
polynomials of degree d up to swapping of variables. 

In degree 2, the homogeneous polynomial (x + y) 2 = x 2 + 2xy + y 2 is sharp. 
This polynomial is the only known case of a sharp polynomial with 3 terms of 
top degree. The only other polynomial (up to swapping of variables) in degree 
2 is the so-called Whitney map x + y(x + y) = x + xy + y 2 . 

We summarize results for the number of polynomials in H(2,d) where d = 
2N - 3, d = 2N - 4, and d = 2N - 5 in Table [U Here we do not identify 
polynomials up to swapping of variables. The first row is the number of sharp 
polynomials for odd degrees, the second row is the number of sharp polynomials 
in even degrees. We note that verifying the first two rows of the table essentially 
amounts to proving Theorem 12.21 For degree d = 2N — 6, with N > 4, there 
always exists a one parameter family of polynomials in "H(2, d) and hence the 
number is always infinite. It is unknown if there are any one parameter families 
when degree is d — 2N — 5, that is, if the number of such polynomials is 
infinite for some N. A negative answer to this question is equivalent to proving 
Proposition 15.31 for all degrees. 

We pose several open questions about the sequences corresponding to rows 
of Table [2] Denote these sequences by ojv, 6jy, and cat corresponding to the 
first, second, and third row respectively. Is bounded? Both 6jy and cn are 
unbounded, but we know little about their rates of growth. Is or 
bounded? Is 6jy or cjv monotone? If Theorem 12.21 is true for all degrees, then 
6jv must be monotone, as we can derive a formula for 6j\r in terms of the first 
row of the table (with the exception of 63). That is, = 2(a2<Z7v_i + a^ajy^i + 
■ ■ ■ + aN-ia2)- None of the three sequences previously appeared in The On-Line 
Encyclopedia of Integer Sequences jSloj . and were entered as A143107, A143108 
and A143109. 

The results above, together with the results given in §[7] and §[8] also suggest 
the following conjectures. 

(i) Up to swapping of variables, fd is the unique sharp polynomial in "H(2, d) 
for infinitely many odd degrees d. 

(ii) All sharp polynomials in even degrees greater than two are constructed 
from sharp polynomials in smaller odd degrees by the procedure in ([3]). 
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Table 1: Complete list of sharp polynomials of odd degree up to 17. Degrees in 
which uniqueness holds are marked with an asterisk. 

Sharp polynomials 



x + y 

x 3 + 3xy + y 3 

x 5 + 5x 3 y + 5xy 2 + y 5 

x 7 + 7x 3 y + 14x 2 y 3 + 7xy 5 + y 7 
x 7 + 7x 3 y + 7x 3 y 3 + 7xy 3 + y 7 
x 7 + \x h y + \xy + \xy° + y 7 

x 9 + 9x 7 y + 27x b y 2 + 30x 3 y 3 + 9xy 4 + y 9 

x 11 + llx 9 y + Ux 7 y 2 + 77x 5 y 3 + 55x 3 y 4 + llxy 5 + y 11 
x 11 + llx 5 y + llx 5 y 5 + 55x 4 y 3 + 55x 3 y 5 + llxy 5 + y 11 

x 13 + 13x n y + 65xV + 156x 7 y 3 + 182x 5 y 4 + 91x 3 y 5 + 13xy 6 + y 13 
x 13 + 13x n y + 65x 9 y 2 + ^-x 7 y 3 + f x 3 y 3 + ^x 3 y 7 + 13xy 6 + y 13 
x 13 + f^y + if x*y 2 + if *V + f b xy + ifxy 6 + §xy" + y 13 
x 13 + f^ n y + if x 9 y 2 + if x 7 y 3 + %xy + if xy 6 + §xy" + y 13 

x 15 + 15x 13 y + 90x n y 2 + 275x 9 y 3 + 450x 7 y 4 + 378x 5 y 5 + U0x 3 y 6 
+ Ibxy 7 + y 15 

x 15 + 140x 9 y 3 + 15x 7 y + 420x 7 y 4 + 15x 7 y 7 + 378x 5 y 5 + 140x 3 y 6 
+ 15xy 7 + y 15 

x 17 + 17x 15 y + 119x 13 y 2 + 442x n y 3 + 935x 9 y 4 + 1122a; 7 ?/ 5 
+ 714x 5 y 6 + 204x 3 y 7 + 17xy 8 + y 17 



Table 2: Number of polynomials in the top 3 degrees for each N. 



degree 
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5 
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10 
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4 
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4 
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4 
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10 


24 


32 


56 
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d = 2N - 
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11 


38 


88 


198 
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? 
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3 Form of sharp polynomials 



In this section we prove new results about the form of sharp polynomials. These 
results are proved by extending the graph theoretic proof of the inequality 
d < 2N — 3 from DKR03 . We therefore sketch the main ideas of that proof 
as we will need them. Our primary application of these results is to reduce 
the computation time of our computer code, but they are also of independent 
interest. The main result of this section is the following lemma. 

Lemma 3.1. Let d be an odd integer and let f £ H(2,d) be sharp. Then, 

f(x,y) = x d + y d + (lower order terms). (4) 

We remark that in general it is not too hard to prove that any / £ %{2,d) 
must have at least two terms of degree d, and such a statement can be generalized 
to higher dimensions as well. Also the reader should notice that the lemma 
cannot possibly hold for d even, see ((3]). In fact, there exists no polynomial in 
H(2, d) of even degree of the form (@|. Writing f(x, y) — 1 = (x + y — l)q(x, y), 
we notice that the top degree terms of / must be divisible by (x + y) and x d + y d 
is only divisible by x + y if d is odd. 

Proposition 3.2. Suppose f £ H(2, d) and f(x, y) = x d + y d + g(x, y), where 
g is of degree less than d. Then g(0,y) = g(x,0) = 0. 

In other words, if / is of the form ((4|), then the lower order terms involve 
only mixed terms. By mixed terms we mean terms involving both x and y. 

Proof. Notice that /(1,0) = 1. Then it is obvious that g(l,0) = 0. As the 
coefficients of g are positive, we get g(x, 0) = 0. □ 

A further corollary of Lemma |3. II also reduces the search space: 

Lemma 3.3. Let d > 1 be an odd integer and let f £ TL(2,d) be sharp. Then 
at least one monomial of degree d — 1 has a nonzero coefficient. 

Proof. Write f(x,y) — 1 = (x + y — l)q(x,y) for some q. Write q = qd-i + 
qd-2 + • • • + <7o as the homogeneous decomposition of q. If / does not have any 
monomials of degree d — 1 then (a; + y)qd-2 — q<i-i = 0. Hence, the top degree 
terms of / must be divisible by (x + y) 2 . By Lemma [57X1 the top degree part of 
a sharp polynomial is x d + y d , which is not divisible by (x + y) 2 . □ 

Further analysis of the proof of [DKR03] reveals the following lemma, which 
further reduces the search space. 

Lemma 3.4. Let d > 1 be an odd integer and let f £ H(2, d) be sharp. Then f 
does not contain terms of the form x^y^^i for even j . 

Before we prove Lemma 13. II and Lemma 13.41 we need to set up the terminol- 
ogy of [DKR03 and restate some of their results. First we have the following 
characterization of homogeneous polynomials in H(2,d). More general related 
results were proved in |Rud84] and |D'A93j . For convenience of the reader we 
prove the following special case, in the setting of that applies to this work. 
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Proposition 3.5. Let h £ H(2,d) be homogeneous. Then h(x,y) = (x + 

Proof. Note that (x + ?/) = h(x,y) when x + y = 1. Any point in the first 
quadrant of ]R 2 can be written as (tx,yx) where x + y = 1 an t > 0. By 
homogeneity, (te + i?/) d = t d (x + y) d = t d h(x, y) = h(tx, ty). Two polynomials 
equal on an open set are equal everywhere. □ 

Multiplying any lower degree part of a polynomial in ~H(2, d) by (x + y) k does 
not get us out of the space. Hence, for / £ H{2, d) we write the homogeneous 
decomposition f = f d + f d _ 1 + •• • + /o- Using Proposition 13 . 51 we get 

(x + y) d = f d (x, y) + (x + y)f d ^ + (x + yff d . 2 + ■ ■ ■ + (x + y) d f Q . (5) 

Therefore, every polynomial in H(2, d) is constructed by starting with (x + y) d , 
partitioning it into two parts and dividing one by (x + y) , then repeating the 
process. This operation is called undoing. 
Write 

f(x,y) - 1 = (x + y- l)q(x,y). (6) 

We study the coefficients of q. In particular we write the Newton diagram for 
q where we ignore the size of each coefficient and only write P, TV or for 



positive, negative or zero respectively. For example, when / 



3xy + y a 



then q — y+y — xy + x +x + l and the diagram (including the corresponding 
monomials) is 



(7) 



We have highlighted the entries corresponding to terms of /. For every entry 
in the diagram we define the 2 by 2 submatrix that includes the entry itself, 
the entry just below, and the entry to the left. If the submatrix is any of the 
following, then we say the entry is a sink. 



X 3 
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P N 
* P 
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* P 



P N 
* 



N 
* 



P 

* P 





* P 



P 

* 



(8) 



Each sink in the diagram for q must correspond to a nonzero positive term in 
/. There may be more positive terms in / than there are sinks, but not the 
other way around. In 0, the sinks are marked in bold. One can make the 
corresponding definition of a source, which would force a negative term. There 
must therefore be at most one (and in fact exactly one) source corresponding 
to the —1 in f{x 7 y) — 1. 

From equation ([5]) we see that the corresponding diagram for any / £ H(2,d) 
is obtained by starting with the diagram for (x + y) d and successively changing 
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P's to iV's or O's. The diagram for (x + y) 3 is 





P 

P P 

P P P 



(9) 



We remark that diagrams we obtain by changing P's to TV's or O's need not 
correspond to polynomials in H(2, d). 

D'Angelo, Kos and Riehl prove that the minimum number of sinks one can 
obtain by this procedure is exactly f^i^] . It turns out that there exist poly- 
nomials in H(2, d) with precisely this many nonzero terms. Thus for a sharp 
polynomial in 7^(2, d) the nonzero terms correspond exactly to sinks in its New- 
ton diagram. We can thus easily prove the following version of Lemma l3~T1 This 
proposition was essentially proved in [DKR03] but not stated explicitly. 

Proposition 3.6 (D'Angelo-Kos- Riehl). Let f £ H{2,d) be sharp. Then for 
some k,m < d, 



Proof. Note that there must be exactly one source, hence the lower right entry 
in the diagram for / must be a P. In the bottom row, we note that to only have 
sinks except for the one source, we must have a row of some number of P's and 
then all zeros. That means that there is exactly one sink on the bottom row. 
Similarly there is exactly one sink in the leftmost column. Since / is sharp, 
these sinks correspond exactly to nonzero terms in /. Thus there is at most one 
pure term in x and at most one pure term in y. By plugging in x = 1, y = 0, 
and vice versa, we get that the coefficient of both terms must be 1. □ 

In the procedure of undoing, we start with the diagram for (x + y) d , which 
has d + 1 sinks along the main diagonal and we change P's to jV's or O's. Of 
course we can change only the entries corresponding to terms of degree d — 1 or 
less as the diagram corresponds to the q in ([6]) , which is of degree d — 1 . We call 
a diagram D' an ancestor of D if D has less than or equal number of nonzero 
entries (P's or N's) than has D'. 

We note what can happen in this procedure to the sinks of the diagram 
corresponding to [x + y) d , which we call Dh- The sinks may move leftward in 
rows, or downward in columns, in which case the number of sinks is unchanged. 
Sinks can also be created, two sinks can coalesce into one or no sinks, or a 
sink can disappear. The main idea of the proof is essentially the following 
result which we state as a lemma, and which is proved in [DKR03] . By a 
procedure of getting to a diagram D of / we mean a finite sequence of diagrams 
Dh = Da — > Di — > ■ ■ ■ — > Dm = D such that Dj is an ancestor of Dk whenever 
j < k. All diagrams are ancestors of D and Dh is the ancestor of all the diagrams 
in the sequence. 

Lemma 3.7 (D'Angelo-Kos- Riehl). Let D be a diagram corresponding to a 
sharp f £ T-L{2 1 d). Then there exists a procedure (as described above) such that 



f(x,y) 



x ' + y m + (mixed terms). 



(10) 
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each diagram in the procedure has a unique source at the origin. Furthermore, 
all coalescence happens on the diagonal corresponding to terms of degree d — 1 . 
In particular, there exists an ancestor diagram D\ with the same number of 
sinks as D, and with sinks corresponding only to terms of degree d — 1 and d. 

Proof. For completeness we sketch one possible proof of this lemma, which is 
slightly different from |DKR03j . The tedious details require checking a finite 
number of cases and are left to the reader. We work in reverse, that is, we start 
with D and work towards D\. Take D and notice that we can change O's to P's 
or TV's without increasing the number of sinks. This can be done in such a way 
that sinks in move up and to the right or disappear. It is not hard to check now 
that by setting certain iV's to P's we can have sinks move up and to the right 
until they reach the diagonal corresponding to degree d — 1 or disappear. This 
diagram is the D\ we seek. □ 

We are now ready to prove Lemma 13.11 Intuitively the idea of the proof is 
that if all the coalescence happens on the diagonal corresponding to terms of 
degree d — 1 then the coalescence of sinks must happen in matched pairs. An 
extra sink is left over on each side of a row of such matched pairs. To minimize 
the number of sinks, there must be at most two such extra left over sinks and 
they must be the ones corresponding to x d and y d . 

Proof of Lemma WA\ Suppose / € %{2, d) is sharp and d is odd. We now know 
that during the procedure of getting the diagram for /, we have passed through 
a diagram D\ with sinks only on the diagonals corresponding to degrees d — 1 
and d, and having the exact same number of sinks as the number of nonzero 
terms in /. 

Let us start with Dh- It is easy to see that except for the first and last entry 
on the degree d — 1 diagonal, maximum coalescence happens when we change 
every other P on the diagonal to N. Changing to does not remove the sinks 
of degree d and hence does not create coalescence. Changing the first or the last 
entry to N does not create any coalescence. Since d is odd, we note that the 
degree d — 1 diagonal of the diagram consists of PNPNP . . . PNP. Any other 
arrangement has too many sinks. 

The D\ diagram has sinks for x d and y d and no other sinks of degree d. We 
can no longer create or lose any sinks on the diagonal of degree d. We know that 
as / is sharp, sinks correspond exactly to terms in /. We know that there are 
at least two terms of degree d. Applying Proposition 13. 61 gives the result. □ 

Proof of Lemma \3J\ We analyze the proof further. We again find the diagram 
D\ that has PNPNP . . . PNP as the degree d — 1 diagonal. If we create a sink 
by changing one of the P's to a zero or an N we note that we have created a 
source, which is not allowed. □ 
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4 Finding sharp polynomials as a linear problem 



We fix d and we treat the coefficients of polynomials as variables. Hence, we 
can treat the space of polynomials of degree d or less as for some large K. 
Suppose 



That means all the non-constant coefficients of h(x) = p(x, 1 — x) must be zero 
and the constant coefficient must be 1. We get a linear system of d equations in 
K variables and one affine equation. If we instead let p(x, 1 — x) = c and let c 
be a variable, then we get a linear system of d+ 1 equations in K + 1 variables. 
Whenever we find a solution p to p(x, 1 — x) = c, c 0, we get a solution to 
p(x, 1 — x) = 1 by rescaling. 

Now that we know that p 6 H(2, d) are solutions of a linear system of equa- 
tions, we prove some useful propositions. We will need the following proposition 
from [Lcb07b] , which we reprove here as the idea of the proof is important in 
the next section. A generalized version of this proposition was given in [DL09b . 

Proposition 4.1. If there exists a continuous map t i— >• pt € T-L[2, d) (a one 
parameter family) and further that N{j> t ) is constant for t in some open interval 
I . then pt is not sharp for t G I . 

Proof, p G 1-1(2, d) are solutions of a linear system, hence if there is a family, 
then there is a straight line with the same property. We pick two polynomials 
ip and ip such that pt := tip + (1 — t)tjj G T-L(2, d) and the number of terms in p t 
is constant some small interval I. We restrict our attention to a closed interval 
J where pt has nonnegative coefficients. It is not hard to see that J must be 
bounded and we could by rescaling assume that J = [0, 1]. Further, we note that 
the number of nonzero coefficients of po and p\ must be smaller than the number 
of coefficients in p t for t G (0, 1). Obviously / C (0, 1). It remains to show that 
Pq or pi are in "H(2, d). They could conceivably be of lower degree than d, but 
they cannot both be such since pt is a convex combination of them. □ 

We define a support as a subset of the set of monomials of degree at most d. 
The support of a polynomial p is the set of monomials with nonzero coefficients. 

Proposition 4.2. No two sharp polynomials in 1-1(2, d) have the same support. 
In particular, there can be at most finitely many sharp polynomials in H(2,d). 

Proof. If p and q have the same support, take the combination tp + (1 — t)q, 
which is a one parameter family of same support. By Proposition 14.11 cither 
p = q or neither p nor q can be sharp. □ 

Corollary 4.3. The coefficients of sharp polynomials are rational. 
We omit the proof as we never use this result. 




(11) 



j+k<d 



The condition that p(x, y) = 1 on x + y is equivalent to 

p(x, 1 — x) = 1. 



(12) 
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5 Linear algebra method 



For j, k satisfying < j ' + k < d — 1, we define polynomials 
b jtk (x, y) := x J y fe - xV(z + y) d ~ j ' k . 



(13) 



Each equals on the line x + y = 1, and they are linearly independent. By 
counting dimensions it is not hard to see that the polynomials (ITBl together 
with (x + y) d span the space of all polynomials that are constant on the line 
x + y = 1. (|13p helps us construct the actual system of linear equations to find 
sharp polynomials in H(2,d) in a relatively simple way. 

We decompose p e H(2, d) using the basis (fTU)) together with (x + 



p(x,y) = (x + y) d + ^2cj, k bj,k(x,y). 

3,k 



(14) 



We note that c J: fc are precisely the coefficients of p and we also note that since 
p{x, y) = 1 on (x + y) = 1, then the coefficient of (x + y) d in the decomposition 
must be 1. Hence the coefficients of p degree d are affine functions of the 
coefficients of degree d — 1 or less. If we treat the coefficient of (x + y) d (and 
thus the value of p on the line x+y = 1) as a variable we find that the coefficients 
of p are linear functions. 

We will always assume that co ; o = even if looking for nonsharp polynomials 
in H{2, d). If p(0, 0) ^ then p(a 1 '^ p( ° ) ,0) € H(2, d). If we are only looking for 
sharp polynomials, then Cq.o will always be zero. 

We construct a matrix A with each column corresponding to one b^k and 
one column corresponding to (x + y) d . Each row represents one monomial of 
degree d. This matrix when applied to the vector (. . . , Cj t k, . . . , 1)* produces 
a vector of the degree d coefficients. For illustration, suppose that d = 3 and 
we order our monomials of degree 2 or less as (x, y, x 2 , xy, y 2 ) and order the 
monomials of degree 3 as (x 3 , x 2 y, xy 2 , y 3 ). Then 



.4 



-10-10 1 
-2 -1 -1-10 3 
-1 -2 -1-13 
0-10 0-11 



(15) 



Now A(cio,coi,c 2 o,cii,co2, 1)' = (030,221,012,003)'. 

We need an algorithm to test if a given support is a support of a sharp 
polynomial in H(2,d). We divide the support into the degree d part and the 
lower degree part. We pick all the monomials of degree d— 1 or less in the support 
and pick out the corresponding columns in the matrix A, plus we always take 
the column corresponding to (x + y) d . Then we pick out the rows corresponding 
to the monomials that do not appear in the support. We get a submatrix A' 
and compute its nullspace. If A' has empty nullspace, no polynomial with such 
a support can vanish on the line x + y = 1 . 

Any vector in the nullspace of A' represents a set of coefficients of degree 
d — 1 or less in a polynomial that vanishes on x + y = 1. We apply these 
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coefficients to A to obtain the degree d coefficients. The coefficient of (x + y) d 
must not be zero if all the coefficients of the polynomial are to be positive, as a 
polynomial that is zero on x + y — 1 must have coefficients of both signs. 

We claim that the nullspace must be of dimension exactly 1 if it corresponds 
to a sharp polynomial in H(2, d). If the dimension is more than one we would 
obtain a family of polynomials with the given support by Proposition 14.11 the 
polynomials cannot be sharp. 

Therefore we compute the nullspace of A'. If it is of any other dimension 
than 1, we are done. If it is of dimension 1, apply the corresponding vector to 
A and test all coefficients for being nonnegative. 

As a simplification we note that except for the last column, A (and hence 
A') consists of nonpositive numbers, and the last column is positive. Hence if 
A' contains a row of the form (0, 0, ... , 0, c) for some constant c, the nullspace 
cannot contain a nonzero vector with only nonnegative entries. 

As an example we take the A for degree 3 as given above. Suppose that we 
wish to test the support xy,x 3 ,y 3 , that is 

(16) 

The nullspace of A' is exactly one dimensional and (3, 1)' spans this space. 
Therefore p(x, y) = (x + y) 3 + 3&n or p(x, y) = 3xy + x 3 + y 3 . Further compu- 
tation shows that 3xy + x 3 + y 3 is the only sharp polynomial in "H(2, 3). 

When searching for sharp polynomials we apply the following simplifications. 

(i) When the degree d is odd, the terms of degree d are precisely the terms 
x d + y d . See Lemma l3~Tl 

(ii) Exactly two pure terms occur. See Proposition 13.61 

(iii) When the degree d is odd, at least one term of degree d — 1 appears. See 
Lemma 13.31 

(iv) When degree d is odd, terms of the form x^y^ 1 ^^ do not appear for even j. 
See Lemma [331 

(v) From any two terms of the form x k y m+1 and x k+1 y m , only one occurs. 
See Remark |5 . 1 1 below . 

Remark 5.1. Suppose our polynomial contains c k ,m+ix k y m+1 + Ck+i, m x k+1 y m , 
for nonzero Ck, m +i and Ck+i, m - Assume without loss of generality that Ck. m +i < 
c k+1 , m . Then c k , m+1 x k y m+1 + c k +i, m x k+1 y m = c k , m+1 (x + y)x k y m + {c k+1 , m - 
c k,m+i)x k+1 y m ■ We can replace x + y with 1, to obtain a new nonequivalent 
sharp polynomial. Once we have found the second one we would have found 
the first as well. No sharp polynomials with this configuration of monomials 
have been found so far. If Proposition 15.31 below is true for all degrees then 
no such sharp polynomials actually exist. Let / be the polynomial containing 
two monomials as in (jvj), and g be the polynomial obtained by dividing out an 
(x + y). We would obtain a contradiction by considering (1 — t)f + tg. 



A' = 



-1 3 
-1 3 
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The algorithm in this section has been implemented using the Genius soft- 
ware version 1.0.2 |Leb07aj . Parts were also implemented in plain C using the 
GMP library [GMP07] . To reduce computations using large integer arithmetic, 
we used modulo arithmetic first to check if the matrix A' is nonsingular. A very 
small prime, p = 19, was sufficient to eliminate vast majority of cases. 

The tests were run on a recent Intel 3GHz CPU. For d = 11 the time used 
was 0.18 seconds, for d = 13 the time was 17 seconds, for d = 15 the time 
was 33 minutes, and finally for d = 17 the time was 77 hours. From these 
timings it appears that the complexity of this method grows faster than the 
mixed linear programming method described in § [5] However, at least with the 
current implementations, the method of this section is faster for small degrees. 

We can also use this method to find even nonsharp polynomials. In this case 
we compute the nullspace, which is now possibly more than one dimensional. It 
is no longer easy to find the subspace generated by nonnegative vectors. We get 
many families of polynomials that have negative coefficients. A simple heuristic 
can eliminate most such families. The rest are easy to sort through by hand. 

In |DL09bj . D'Angelo and the first author proved the following theorem. By 
a fc-dimensional family of polynomials we mean a fc-dimensional polytope in the 
parameter space. Recall that N(f) is the number of distinct monomials of /. 

Theorem 5.2 (D Angelo-Lebl). Let T C W(2,d) be a k- dimensional family, 
then for any f € J- 

d<2(jV(/)-*)-3. (17) 

The bound (|17p is not sharp. By using the algorithm to find all polynomials 
in H(2,d) as above, we get a computer assisted proof of the following improve- 
ment in a special case, and this result is sharp. That is, no better inequality is 
possible for d < 9 and a 1-dimensional family. 

Proposition 5.3. Let J- C H(2,d) be a 1-dimensional family, d < 9, then for 
any f € T 

d<2N(f)-6. (18) 

6 Mixed linear programming method 

Another approach to computing sharp polynomials in "H(2, d) involves constraint 
satisfaction of mixed-integer programs. That is to say, some variables will take 
on integer (actually 0-1) values and others will be continuous. We will ap- 
ply classical methods, using branching and a naive form of cutting planes for 
handling the integer variables. Good background references are |Dan63j and 
|Sch86j . 

We now describe in brief the setup of such problems. We are given an odd 
degree d. We prescribe that the minimal number of nonzero terms is i.e. 
the polynomial is sharp. From theory just presented we know several constraints 
on such polynomials. Therefore we also prescribe all the constraints listed in 
the last section. 
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We now assign two variables to each of the monomials not ruled out. One 
variable will record whether that monomial goes into the polynomial; it takes on 
values or 1. The other is constrained to be nonnegative and corresponds to the 
coefficient of the monomial in the sharp polynomial we attempt to construct. 
These latter variables satisfy equations that arise from the identity p(x, 1 — x) — 
1 = 0; we set the coefficients of terms in every degree to zero. 

As shown earlier, one can multiply terms of the sharp polynomial with pow- 
ers of (x + y) to obtain the homogeneous polynomial in H(2, d) of degree d (and 
indeed, we get the same equations as from zeroing coefficients of p(x, 1 — x) — 1). 
This procedure in fact gives upper bounds on each coefficient: if x^y k has co- 
efficient Cj ; fc then one can easily show that < Cj t k < min{(^), (^)}. With 
a bit more work one can deduce a generally stronger inequality: for each 
Q < m < d — j — k we have Cj t k < im) / ( d ~m~ k ) ■ Nonnegativity is imposed 
by the requirements of our polynomials. The second inequality is of significant 
interest, because it allows us to relate the continuous variables to their discrete 
counterparts. Specifically, if we call the corresponding 0-1 variable bj^ and call 
the minimum of the binomial quotients m,j } k, then we have Cj t k < m j,kbj } k- 

To simplify the computations we may impose a few other restrictions; gen- 
erally speaking, the more inequality restrictions we impose, the faster the com- 
putation runs. As before, we can insist that there be no pair of neighbors of 
equal degree in the polynomial. That is, we can impose that the sharp polyno- 
mial does not contain Ck, m +i x y m + Ck+i, m x k+1 y m , for nonzero Ck, m +i and 
Cfc+i, m - See Remark 1 5. II 

Another restriction is to insist that the sum of coefficients from the left side 
of the Newton diagram (that is, monomials with deg(x) > deg(y)) be larger 
than or equal to the sum from the right side. This asymmetry is valid since we 
do not care about equivalent polynomials obtained by exchanging variables. 

The actual code begins by setting up linear equations and inequalities based 
on the discussion above. We solve the equations, thus eliminating some vari- 
ables. We use this result to adjust the inequalities accordingly. The implemen- 
tation is a standard branch-and-cut mixed linear programming code. We branch 
on the variables constrained to be integral. 

The Mathematica code to do what we have described in this section occupies 
about 60 lines. It handles the case d = 9 in about 1.3 seconds, d = 11 in 24 
seconds, d — 13 in 9 minutes, d = 15 in 4.6 hours, and d = 17 in 186 hours. 
These timings are on a fairly recent CPU operating at 3.2 GHz, running version 
6.0.2 of Mathematica, with settings to use the COIN-LP library }LH03] to solve 
the relaxed linear programs. They are of course also dependent on the extent 
to which the authors have found algorithm simplifications based on the theory 
presented for these polynomials. A better understanding of terms that must or 
must not arise in such polynomials would almost certainly lead to algorithmic 
improvements . 
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7 Uniqueness of sharp polynomials 



Partial information about when uniqueness fails can be summarized in the the 
following theorem. This theorem combines the results of this paper with the 
results of [DL09aj . The fact that uniqueness holds when d — 9 was given without 
proof in [DL09a in anticipation of the present paper. 

Theorem 7.1. Uniqueness holds when d = 1, d = 3, d = 5, d = 9, and d = 17. 

Uniqueness fails in the following cases: 

(i) Suppose d is even. Then uniqueness fails for all d. 

(ii) Suppose d is congruent to 3 mod 4. Then uniqueness holds for d = 3 and 
fails for d > 7. 

(Hi) Let k be a positive integer. Uniqueness fails for d of the form 

_ (7 + 4V3) fc + (7 - 4V3) fc 
2 

I.e. d = 7, 97, 1351, 18817, 262087, . . . 

(iv) Suppose d > 1 and d is congruent to 1 mod 6. Then uniqueness fails. 

The first case not handled by this theorem is d = 21. It is computation- 
ally infeasible to completely test this case with the algorithms we have so far. 
Nonuniqueness in the theorem is proved by an explicit construction of new sharp 
polynomials that is sketched out in the next section. We can test these construc- 
tions more generally using a computer and get further results on degrees where 
uniqueness fails. We get a computer generated proof of the following proposi- 
tion, which contains more information than Theorem 17. II for small degrees. 

Proposition 7.2. Uniqueness fails for all degrees d < 149 not contained in the 
following list: 

1, 3, 5, 9, 17, 21, 33, 41, 45, 53, 69, 77, 81, 93, 105, 113, 117, 125, 129, 141, 149 (19) 

More precisely, (| 19[) lists all degrees d < 149 where the procedure of § [S| fails to 
produce a sharp polynomial besides the group invariant one. 

We have run the computer code for degrees up to 2 9 + 1 and the sequence 
above does not appear to thin out very rapidly. It seems reasonable to conjec- 
ture that the sequence is infinite. Since we also know that, at least up to degree 
17, that the construction of § [8] gives all sharp examples, it is also reasonable 
to conjecture that the sequence of degrees for which uniqueness holds is infi- 
nite. The On-Line Encyclopedia of Integer Sequences Slo j did not include the 
sequence in (fl~9"|) . It is now entered as A143105. The partial known sequence 
of degrees where uniqueness holds was entered as A143106. The encyclopedia 
does not contain any sequence that is a subsequence of (fT9|) and starts with 
1,3,5,9,17. 
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Note that the beginning of the list (up to 33, with the exceptions of 1 and 21) 
are degrees of the form 2 k + 1. It would be incorrect to assume that uniqueness 
holds for such degrees. In fact, uniqueness fails for example in degrees 2 6 + 1, 
2 8 + 1 and 2 9 + 1. 

8 Construction of new sharp polynomials 

Except for the even degree case, all the new noninvariant sharp polynomials con- 
structed for the proof of Theorem 1 7 . 1 1 arise in a similar way. For the nonuniquc- 
ness in the even degree case see ©. There we construct a sharp polynomial 
of even degree d = d\ + d 2 from two sharp polynomials of odd degrees d\ and 
g?2- We observe that the number of possible sharp polynomials goes to infinity 
for even degrees. We also remark that using this method we can also construct 
group invariant sharp polynomials by taking d\ = d-z = i and using a group 
invariant polynomial of degree ^ in ((3)) . 

Suppose that d is odd. Take some even m < d and look at f m . Write 

f m (x,y) = f m (x,y)-y m , (20) 

Now notice that f m = 1 on x + y = 1 , hence on x + y = 1 we have 

f m (x,y) = l + y m . (21) 

If we can find a constant c and a monomial x^y k , such that cx^y k f m (x,y) has 
at least two terms common with fd(x,y), we can "replace" cx^y k f m (x,y) with 
cx : 'y k (l + y m ). More explicitly, we can write 

f(x,y) :=f d (x,y)+cx^y k (l + y m -f m (x,y)) 
= fd(x,y) - cx 3 y k {f m {x,y) - 1). 

From the right hand side, it follows that f(x,y) = 1 on x + y = 1. If we can 
show that / has all positive coefficients, then we are finished. 

As an example, take d — 7. Then fj = x 7 + 7x 5 y + 14x 3 y 2 + 7xy 3 + y 7 and 
f2 = x 2 + 2y — y 2 . Then we get the new sharp polynomial 

f(x, y) ■= h(x, y) - 7x 3 y(f 2 {x, y) - 1) 

= x 7 + 7x b y + Ux 3 y 2 + 7xy 3 + y 7 - (7x 5 y + 14x 3 y 2 - 7x 3 y 3 - 7x 3 y) 
= x 7 + 7x 3 y 3 + 7x 3 y + 7xy 3 + y 7 . 

(23) 

All sharp polynomials known to the authors are obtained by (|22j) in one 
or more steps. In particular, we know all noninvariant sharp polynomials of 
odd degree d < 17 are constructed this way. When d = 13 then two steps 
are required, that is we need to repeat the procedure (f22j) twice, to obtain one 
of the noninvariant sharp polynomials. To find new polynomials this way, we 
look for certain ratios of terms in the possible coefficients of fd, and then verify 
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that the procedure does not introduce negative terms. Each known noninvariant 
sharp polynomial gives rise to an infinite sequence of degrees satisfying a certain 
Pell equation (see [DL09a ) where the same construction applies. It is still 
necessary, however, to check that all the coefficients in the new polynomials are 
nonnegative. If all the coefficients are nonnegative, we get an infinite sequence 
of new sharp polynomials in different degrees. For a few specific cases this 
is essentially what was done in [ DL09aj to prove the nonuniqueness parts of 
Theorem 17.11 The computations, while elementary, quickly become long and 
tedious. The Pell equation is degenerate in one very specific case when d = 1 
(mod 6). Otherwise, the sequence obtained is very sparse and "thins out" very 
quickly as the degree rises. See |DL 09a for more information. 

We have formulas for the coefficients of terms of fd, that is, we know that 
except for s = 0, the coefficient of x d ~ 2s y s in fd is f ^ee ^ or exam_ 

ple |D'A93j . Hence, it is not hard to check via computer for degrees where 
the construction produces new examples. We have run simple computer code 
to try this procedure for all possible parameters for d < 513. A partial list of 
degrees where the procedure fails to produce new sharp polynomials is given in 
Proposition 17.21 

9 Background in complex geometry 

As stated earlier, the primary motivation for this work comes from CR geometry. 
Let us describe this connection in detail. For more information about complex 
analysis and CR geometry see the book |D'A93] . There has been much interest 
recently in the CR geometry community in studying the complexity of CR maps 
between manifolds. An interesting model case comes about by studying proper 
holomorphic maps between unit balls in C". Due to the symmetries of the unit 
ball, this problem has nontrivial connections to many areas of mathematics 
including number theory, combinatorics and real algebraic geometry. 

To be more precise, let B„ c C™ be the unit ball and suppose that F : B„ — > 
Bjv is a proper holomorphic map. When F extends to a continuous map of the 
closed ball B„ then the map is proper if it maps the boundary of B„ to the 
boundary of Bjv- It is not hard to see using elementary complex analysis that 
if the target dimension TV is smaller than the source dimension n, then there 
are no proper maps. Alexander |Ale77j proved that when n = N and n > 2, 
then all proper maps are automorphisms of the unit ball, i.e. linear fractional. 
See also the survey |For93] or the book |D'A93j for more information about the 
problem. We will discuss only n = 2 from now on. 

When F extends to the closure of the ball its restriction to the boundary 
defines a CR map from the unit sphere in C 2 to the unit sphere in C N . When 
F is sufficiently smooth up to the boundary then Forstneric |For89] proved 
that F is rational and that the degree of F is bounded in terms of 2 and N. 
Faran |Far82j classified all maps when N = 3. 

D'Angelo has made a systematic study of rational proper maps (see for 
example [D'A030 and classified the polynomial maps |D'A88j in the following 
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sense. If we allow the target dimension to be large enough, all polynomial 
proper maps are obtained by a finite number of operations from the unique 
homogeneous map. D'Angelo conjectured that 

degF<27V-3. (24) 

Further, we will assume that F is a monomial map, that is, every component 
of F is a single monomial. We now change notation slightly. To say that F 
extends to the boundary is to say that ||F(z)|| 2 = 1 whenever ||z|| 2 = 1, where 
||-|| is the standard euclidean norm on and C 2 respectively. When F is a 
monomial map, then we can replace |zi| 2 by the real variable x and \z2\ 2 by the 
real variable y. ||-F(z)|| then becomes a real polynomial in x and y of same 
degree as F and with N nonnegative coefficients. Similarly ||z|| becomes x + y. 
Recall that the set of polynomials of degree d with nonnegative coefficients such 
that p(x,y) = 1 whenever x + y = 1 is denoted by 'H(2,d). The inequality (J]) 
proved in [DKR03] therefore proves the conjecture (|24| in the special case when 
the map is a monomial map. Furthermore, when d is odd, the polynomials (p]) 
induce monomial proper maps of balls such that d — 2N — 3. Hence if the 
conjecture (|24| is true for all CR maps, then it is sharp in the sense that the 
bound is the best possible. The maps induced by ([2]) are group invariant and 
hence induce maps of lens spaces. 

A complete classification of the monomial maps is an important first step 
in classification of all CR maps between spheres, and is the main motivation 
for doing the computations in this paper. For example, Faran's result |Far82j 
on the classification of maps from B2 to B3 says that the map is equivalent to 
one of four possible monomial maps. In our language of polynomials in "H(2, d), 
these are the polynomials x + y, x + xy + y 2 , x 2 + 2xy + x 2 and x 3 + 3xy + y . 

The sharp maps arising from fd are group invariant under the action of a 
finite subgroup of the unitary group U(2) generated by § e l for a primitive 
root of unity e. An natural question to ask is to find all other invariant maps. 

Rudin Rud84j proved that a homogeneous proper map of balls is equivalent 
up to unitary transformations to the identity map tensored with itself d times. 
This map is invariant under the action of a cyclic group generated by the matrix 
el, where e is a primitive root of unity and I is the identity. Forstneric |For86j 
observed that only fixed-point-free matrix groups could arise in the context of 
proper holomorphic maps between balls that are smooth up to the boundary. 
Based on this work the second author together with D'Angelo, [Lic92] and 
|DL92j . proved that (up to linear transformations) the only groups T for which 
there exists a T- invariant rational proper map of B2 to B^v are the cyclic fixed- 
point-free matrix groups generated by [ q e m 1 where e is a primitive d th root of 1 
for some odd d, m = 1, 2, or 2m — 1 modulo d. Therefore, the maps generated 
by fd and the homogeneous maps are the only group invariant maps of spheres. 
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